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SUMMARY 


Single-event  upset  of  digital  circuits  due  to  an  ionizing  particle  is  a 
threat  to  an  increasing  portion  of  the  Army’s  electronics.  Extensive  work  in 
modeling  this  effect  has  been  done  at  various  laboratories.  These  modeling 
efforts  have  been  focused  on  field  funnel ing,  and  all  have  neglected  avalanche 
multiplication,  even  though  one  experimental  paper  from  the  Naval  Research 
Laboratory  (NRL)  has  shown  the  importance  of  avalanching.  This  report  veri¬ 
fies  that  importance  by  obtaining  good  agreement  between  our  computer  cal¬ 
culations  and  NRL  measurements.  Avalanching  may  lead  to  burnout  in  very  high 
speed  integrated  circuits  (VHSIC).  The  present  calculations  are  one  dimen¬ 
sional  and  need  to  be  verified  by  two-dimensional  calculations.  This  report 
is  an  expansion  of  another  work.* 


*Alford  L.  Ward,  Avalanching  In  Single-Event-Upset  Charge  Collection  in  Seed  conductor 
Diodes,  IEEE  Trans.  Nucl.  Scl.,  NS-33  (December  1986),  pp  1552-1559. 
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1.  INTRODUCTION 


When  a  highly  ionizing  particle  produces  a  track  of  dense  ionization  in  a 
semiconductor  and  affects  device  operation  through  charge  collection,  a 
single-event  upset  is  said  to  occur.  There  is  now  extensive  literature 
[1]*  on  this  topic. 

Campbell  et  al  [2]  have  published  measurements  showing  that  avalanche  mul¬ 
tiplication  is  important  in  charge  collection  in  single-event  upsets.  They 
noted  evidence  of  avalanche  multiplication  above  3  V  in  a  silicon  diode  which 
showed  a  breakdown  voltage  of  about  7  V.  Many  papers  (e.g.,  H.  L.  Grubin 
et  al  [3])  have  reported  numerical  studies  of  charge  collection  from 
ionized-particle  tracks,  but  none,  to  our-  knowledge,  have  included  avalanche 
multiplication.  The  emphasis  in  most  of  these  papers  has  been  on  the  field 
Tunneling  effect  [4],  a  result  of  the  extension  of  the  high  field  region, 
which  causes  charge  collection  from  beyond  the  normal  depletion  region  of  the 
device. 

This  paper  presents  our  approach  to  the  avalanche  calculations,  gives 
results  of  calculations  of  avalanche  multiplication  in  silicon  and  gallium 
arsenide  diodes,  and  concludes  with  further  discussion  and  conclusions. 

2.  APPROACH 

Charge  collection  of  tracks  from  an  ionizing  particle  has  been  calculated 
with  the  Harry  Diamond  Laboratories  one-dimensional  (ID)  computer  program 
DIODE.  The  computer  program  has  been  described  previously  [5].  Briefly,  the 
continuity,  current  density  formulation,  and  Poisson's  equations  are  solved 
explicitly  in  their  difference  form.  For  these  calculations,  a  constant  volt¬ 
age  was  applied  through  a  series  resistance  of  50  a  to  the  diode,  which  was 
shunted  by  a  1-pF  capacitance.  An  area  of  1  *  10"®  cm2  was  chosen  for  most 
calculations  to  determine  the  interaction  of  the  diode  with  the  external  cir¬ 
cuit.  Ideal  Schottky  barrier  diodes  were  simulated  for  all  calculations  ex¬ 
cept  where  a  comparison  is  made  with  a  PN  Junction  diode. 

Uniform  and  equal  densities  of  electrons  and  holes  over  the  track  length 
were  added  to  the  equilibrium  reverse-bias  distributions  as  initial  conditions 
for  the  calculations.  Therefore,  only  the  decreasing  current  portion  of  the 
pulse  is  calculated. 

3.  RESULTS 

3.1  Silicon 

Most  of  the  material  parameters  for  the  silicon  calculations  were 
those  given  by  Sze  [6].  However,  the  important  avalanche  coefficients  were 
taken  from  Van  Overstraeten  and  DeMan  [7].  The  first  calculations  were  made 
for  a  20-jjm  N-type  diode  doped  to  1  »  lO1^  cm"3.  Ideal  abrupt  (Schottky  bar¬ 
rier)  Junctions  were  used  for  most  calculations.  The  calculated  breakdown 
voltage  is  250  V,  and  the  depletion  width  is  19  ym  at  breakdown.  A  negative 
differential  resistance  (second  breakdown)  is  calculated  at  327  V. 


'References  appear  on  p  22. 


Calculations  were  made  for  an  ionized  track  density  of  1  *  1011*  cm“3 
in  order  to  eliminate  appreciable  space-charge  distortion.  The  resulting  cur¬ 
rent  density,  with  and  without  avalanche  multiplication,  is  plotted  in 
figure  1  for  an  applied  voltage  of  200  V  and  a  20-ym  track  length.  For  fully 
saturated  velocities  of  1  *  10?  cm/s  and  no  diffusion,  the  current  density 
would  drop  linearly  from  320  A/cm2  to  zero  at  0.2  ns.  Actually,  the  field 
drops  to  near  zero  at  17.4  ym,  and  the  hole  velocity  is  not  saturated  over 
about  half  the  diode  width.  The  inset  of  figure  1  shows  the  electron  and  hole 
densities  across  the  diode  at  selected  times  for  the  calculation  without 
avalanching.  It  is  seen  that  the  hole  "front"  is  less  steep  than  the  electron 
front,  because  of  the  higher  hole  velocity  as  the  field  increases  to  the 
left.  Avalanche  multiplication  was  noted  down  to  120  V,  slightly  less  than 
one-half  the  avalanche  breakdown  voltage. 


At  voltages  well  above  the  avalanche  breakdown  voltage,  avalanche 
[8]  (or  IMPATT)  oscillations  are  noted.  These  oscillations  are  comparatively 
large  for  a  track  density  of  2  *  1 0 1 5  cm“3.  Figure  2  shows  the  temporal 
growth  or  decay  for  track  lengths  of  10  and  20  ym.  It  is  seen  that  the 
shorter  track  length,  with  only  one-half  the  initial  charge  of  the  longer 


Figure  1 .  Calculated  current  den¬ 
sity  as  a  function  of  time  for  a 
track  density  of  1  *  10^  **  cm“3  and 
a  voltage  of  200  V.  Dashed  curve 
includes  avalanching,  whereas  dot- 
dashed  curve  does  not.  Inset  shows 
electron  (solid  curves)  and  hole 
(dashed  curves)  density  distri¬ 
butions  for  nonavalanching  case  at 
selected  times  in  picoseconds. 


TIME  (ns) 


Figure  2.  Computed  avalanche  oscil¬ 
lations  are  shown  for  a  20-um  diode 
near  the  second  breakdown  voltage. 
Solid  and  dashed  lines  are  for 
track  lengths  of  20  and  10  pm, 
respectively.  Initial  track  den¬ 
sity  is  2  x  to1 5  cm“3.  Param¬ 
eter  is  applied  voltage. 


track,  has  the  lower  second- breakdown  voltage  (current  increasing  indef¬ 
initely).  This  results  from  a  more  efficient  avalanche  oscillation;  normally 
the  second  breakdown  voltage  decreases  as  the  track  density,  as  well  as  the 
total  charge,  increases. 

The  minimum  voltage  for  avalanche  multiplication,  Vmin,  and  the  second 
breakdown  voltage,  Vg2*  are  plotted  as  a  function  of  track  density  in 
figure  3.  Also  plotted  in  this  figure  is  the  avalanche  breakdown  voltage, 
Vg0,  for  constant-doping,  one-sided,  N-type  diodes  as  a  function  of  that 
doping  density.  The  latter  calculations  pertain  to  diodes  not  punched  through 
at  breakdown;  diodes  which  punch  through  before  breakdown  have  lower  breakdown 
voltages.  The  criterion  for  Vm^n  was  set  at  an  increase  of  collected  charge 
of  about  5  to  10  percent;  a  more  precise  criterion  would  not  be  advantageous 
at  this  time.  These  calculations  were  made  with  various  diode  widths  and 
track  lengths  equal  to  the  diode  width.  For  track  densities  of  lO1^  cm“3  and 
less,  the  diode  width  was  20  pm;  for  track  densities  of  2  *  101?  cm“3  and 
greater,  the  diode  width  was  1  pm.  Intermediate  diode  widths  were  used,  and 
overlapping  track  densities  assured  that  there  was  little  voltage  variation 
with  diode  width  for  non-punched-through  diodes.  This  figure  is  discussed 
further  in  section  4. 


Current  density  as  a  function  of  time  is  plotted  in  figure  4  for  a 
0.5-um  P-type  diode  doped  to  1  *  10^  cm-3  an  initial  track  density  of 
1  *  1018  cm“3t  and  various  applied  voltages.  With  neither  space- charge  ef¬ 
fects  nor  diffusion,  and  assuming  saturated  velocities,  the  charge  would  be 
collected  in  5  ps.  At  just  above  the  minimum  voltage  for  avalanche  multi¬ 
plication,  Vml  -  5  V,  the  charge  collection  time  is  32  ps.  The  charge  col¬ 
lection  time1”  first  decreases  as  the  voltage  increases  and  then  increases  as 
avalanche  multiplication  increases.  The  charge  multiplier,  M,  the  ratio  of 
the  total  charge  collected  to  the  initial  charge  in  the  track,  is  shown  in  the 
inset  of  figure  4.  The  inset  curve  has  the  same  shape  as  the  measured  curve 
of  Campbell  et  al  [2],  namely,  approximately  an  exponential  increase  from  a 
constant  value  at  low  voltages.  The  measured  curves  were  for  similar  doping, 
geometry,  and  track  parameters  as  the  calculations.  Also  plotted  in  the  inset 
of  figure  4  are  the  results  for  calculations  for  a  1-pm  diode. 
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Figure  4.  Current  density  as  a  function  of  time  for 
a  0.5-pm  P-type  diode  doped  to  1  *  1018  cm"3.  Ini¬ 
tial  track  density  is  1  *  1018  cm"3.  Parameter  is 
applied  voltage.  Inset  shows  multiple  of  initial 
charge  as  a  function  of  applied  voltage  for  both 
the  0.5-iim  diode  and  a  1-pm  diode. 


The  inset  curve  data  of  figure  4  have  been  replotted  in  figure  5, 
where  M  -  1  is  plotted  on  a  logarithmic  scale  against  the  voltage.  Also  plot¬ 
ted  in  this  figure  are  the  results  for  a  0.2-pm  diode  with  an  initial  track 
density  of  2  *  1018  cm~3  and  for  a 
0.1 -pm  diode  with  a  density  of 
5  *  1018  cm“3.  in  order  to  plot 
the  experimental  data  given  in  fig¬ 
ure  2  of  Campbell  et  al  [2],  in  our 
figure  5,  it  is  assumed  that  their 
charge  collection  efficiency  (our 
multiplier  M)  was  equal  to  1  at  zero 
applied  voltage.  Their  measured 
data  are  plotted  in  figure  5,  with 
the  further  assumption  that  the 
built-in  potential  is  1  V.  Last, 
calculated  data  for  a  track  density 
of  1  x  io18  cm”  3  and  a  diode  width 
of  1  pm  are  also  shown  in  figure  5 
as  the  dashed  curve.  The  point  for 
V  -  12  V  is  estimated  since  mathe¬ 
matical  calculation  instability  pre¬ 
vented  a  complete  run.  This  is 
addressed  in  section  4.  The  meas¬ 
ured  data  of  Campbell  et  al  [2]  were 
taken  with  a  diode  with  an  epitaxial 
layer  that  was  originally  3  pm,  but 
they  estimate  that  the  layer  was 
reduced  to  1  pm  by  diffusion  during 
the  diode  processing. 


one  is  plotted  on  a  logarithmic 
scale  against  applied  voltage. 
Parameter  for  calculated  curves  is 
initial  track  density.  Diode  widths 
are  given  in  text.  Measured  data 
are  from  Campbell  et  al  [2]. 

Points  at  top  of  figure  are  second 
breakdown  voltages. 


The  inset  of  figure  1  shows 
the  electron  and  hole  distributions 
across  the  diode  during  charge  col¬ 
lection  for  small  (mobile)  space 
charge  conditions.  Figure  6  shows 
the  calculated  distriDutions  of  electrons  and  holes  during  charge  collection 
in  the  P-type  diode  doped  to  1  «  1015  cm"3.  The  initial  track  density  is 
1  x  1018  cm~3  and  the  applied  voltage  is  3  V,  about  the  avalanche  threshold. 
Note  the  steepness  of  the  concentration  gradients  as  electrons  move  to  the 
right  and  the  holes  to  the  left.  The  steepness  results  from  the  higher  field 
and  thereby  higher  velocity  seen  by  the  trailing  carriers.  The  charge  density 
"shoulders,"  especially  noticeable  from  1  to  4  ps  (at  densities  from 
1  x  io11*  to  3  x  1018  cm~3  are  the  result  of  avalanche  multiplication.  The 

fluctuation  of  carriers  about  the  1  x  io18  cm”3  concentration  (shown  only  for 
two  curves)  is  not  physical,  but  a  result  of  incipient  calculation  insta¬ 
bility.  Increasing  the  number  of  calculation  grid  points  removes  the  insta¬ 
bility,  which  is  discussed  later. 


The  field  distributions  corresponding  to  the  distributions  in  figure  6 
are  shown  in  figure  7.  A  logarithmic  scale  is  used  to  show  more  detail.  The 
high  fields  in  the  depleted  regions  are  contrasted  to  the  low  fields  in  the 
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Figure  6.  Carrier  distributions  at 
selected  times  are  shown  for  a  P-type 
diode  doped  to  1  x  10'5  cm“3.  Elec¬ 
tron  densities  are  shown  by  solid  lines 
and  hole  distributions  by  dashed  lines. 
Parameter  is  time  in  picoseconds. 
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Figure  7.  Electric  field  distribu¬ 
tion  corresponding  to  carrier 
distributions  of  figure  6.  Param¬ 
eter  is  time  in  picoseconds. 


the  plasma  region.  The  maximum  field  is  6  *  10^  V/cm  at  0.6  ps.  This  maximum 
field  is  for  the  avalanche  threshold  voltage.  Higher  fields  are  calculated  at 
higher  voltages.  The  time  to  attain  the  maximum  field  assumes  zero  time  for 
track  formation. 

The  current  density  as  a  function  of  time,  corresponding  to  figures  6 
and  7,  wa3  shown  in  figure  *4  as  a  linear  plot  (the  curve  labeled  3  V).  The 
same  data  are  shown  in  figure  8  as  a  log-log  plot  (the  curve  labeled  1  pm). 
Also  shown  in  this  figure  are  results  of  calculations  for  diode  widths  of  from 
0.5  to  5  pm.  Calculations  became  unstable  for  the  wider  diodes  at  times  from 
1  to  10  ps.  The  temporal  current  decay  for  the  same  diodes  is  shown  in  figure 
9  for  an  applied  voltage  of  11  V.  Note  that  the  curves  for  0.5  and  1.0  pm 
cross.  The  other  curves  would  also  cross  if  the  calculations  had  not  become 
unstable. 
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Figure  8.  Current  density  plotted  as 
a  function  of  time  on  a  log-log  scale. 
Track  density  is  1  x  101®  cm “3,  ap¬ 
plied  voltage  is  3  V,  and  diode  width 
in  micrometers  is  given  as  parameter. 
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The  maximum  field  (usually  at  the  cathode,  but  sometimes  at  the 
anode)  is  plotted  as  a  function  of  time  in  figure  10.  Again,  a  log-log  plot 
is  chosen  and  the  data  are  for  the  0.5-ym-width  diode.  The  effect  of  avalanche 
multiplication  is  in  evidence  (at  the  highest  applied  voltages).  It  is  also 
seen  that  the  peak  of  the  maximum  voltage  occurs  at  later  times  for  the  higher 
voltages. 

The  maximum  electric  field  as  a  function  of  time  for  various  track 
densities  is  shown  in  figure  11.  The  diode  widths  for  the  diodes  in  order  of 
decreasing  track  density  were  1.0,  3.0,  15,  and  20  ym.  The  applied  voltages 
were  at  the  avalanche  threshold  and  may  be  obtained  from  figure  3.  The  doping 
density  for  all  diodes  was  1  x  lO1^  cm-3. 


The  maximum  field  across  the  diode  at  any  time  is  plotted  in  figure 
12(a)  as  a  function  of  diode  width.  The  track  density  is  1  x  101®  cm"3f  and 
the  diode  is  doped  to  1  x  1015  cm“3.  For  the  higher  applied  voltages,  calcu¬ 
lations  were  made  with  and  without  avalanche  multiplication.  In  these  cases, 
higher  fields  occurred  with  avalanching  than  without.  At  lower  track  den¬ 
sities,  avalanching  reduced  the  maximum  fields  calculated.  Figure  12(b)  shows 
the  maximum  field  at  zero  (and  at  infinite)  times.  Diodes  below  the  dashed 
lines  are  punched  through  at  the  indicated  voltage;  that  is,  the  field  at  the 
anode  is  zero. 


Figure  10.  Maximum  field  as  a  function 
of  time  for  a  0.5-um  diode  and  a  track 
density  of  1  «  101®  cm" 3.  parameter 
is  the  applied  voltage. 


As  indicated  earlier,  moat 
calculations  have  been  made  simu¬ 
lating  ideal,  abrupt,  Schottky  bar¬ 
rier  diodes.  Calculations  were  also 
made  to  simulate  a  PN  junction 
diode:  the  width  of  the  P  region  was 
0.1  pm  and  the  N  doping  was  1  *  lO1^ 
for  the  1-pm  diode.  The  initial 
track  density  was  1  *  101?  cm“3. 
The  comparison  of  the  two  diodes 
with  an  applied  voltage  of  14  V 
showed  that  the  current  agreed  to 
within  2  percent  at  all  times. 
However  the  maximum  calculated  field 
was  less  for  the  junction  diode,  as 
shown  in  figure  13.  where  the  field 
distributions  are  shown  at  selected 
times  for  the  two  diodes.  The  ap¬ 
plied  voltage  is  Just  above  the  ava¬ 
lanche  threshold,  and  the  charge 
multiplication  is  nearly  equal  for 
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Figure  11.  Maximum  electric  field  as 
a  function  of  time  for  various  track 
densities  (cm"3). 


(a) 


Figure  12.  Maximum  field  (a)  at  any 
time  given  as  a  function  of  diode 
width  and  (b)  at  zero  time.  Diodes 
below  the  dashed  line  are  punched 
through . 
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Figure  13*  Comparison  of  field  distributions  for  (a)  a  Schottky  diode  and 
(b)  a  PN  junction  diode.  N-type  diode  is  doped  to  1  x  io1^  cm~3,  applied 
voltage  is  14  V,  and  initial  track  density  is  1  x  101?  cm"3.  Parameter  is 
time  in  picoseconds. 


both  diodes  because,  although  the  maximum  field  is  less,  the  high  field  region 
is  wider  for  the  Junction  diode.  Maximum  fields  are  attained  at  4  ps  for  the 
Schottky  diode  and  at  5  ps  for  the  Junction  diode. 


Calculated  current  densities  are  plotted  against  time  on  a  log-log 
scale  in  figure  I1!  for  three  initial  track  densities:  1  *  io1^,  i  «  io1?,  and 
1  x  101 8  cm”3.  These  diodes  are  doped  to  1  x  10*5  cm“3  and  the  calculations 
are  made  at  the  avalanche  threshold  voltage,  but  with  the  avalanche  coeffi¬ 
cients  equal  to  zero.  Further  calculations  have  been  made  for  two  and  five 
times  these  track  densities  and  for  at  least  two  diode  widths.  For  clarity, 
two  diode  widths  are  shown  only  for  the  1  *  lO1^  cm"3  track  density.  The 
curves  essentially  coincide  for  the  two  diode  widths  until  the  charge  is  near¬ 
ly  completely  removed  in  the  shorter  diode,  when  its  current  drops  abruptly. 
The  saturated- velocity  transit  time,  Tq,  for  the  8-pm  diode  is  80  ps,  whereas 
the  sharp-drop  transition  time,  Tt,  is  130  ps.  Likewise,  Tq  -  150  ps  and 
Tt  -  700  ps  for  the  15-ym  diode.  When  the  time  ratio  Tt/T0  is  plotted  against 
diode  width  on  a  log-log  plot,  a  straight  line  with  a  slope  of  2  results  for 
each  track  density;  i.e.,  Tt  «  TQ(d)d.  This  is  shown  in  figure  15(a).  Since 
Tq  Is  proportional  to  d,  the  pulse  collection  time  is  proportional  to  d  .  In 
the  same  manner,  a  current  density  ratio  may  be  defined  by  the  ratio  of  the 


current  at  t 


0  to  the  current  at  T 


f 


This  current  ratio  increases  quad- 


ratically  with  the  diode  width  at  each  track  density  calculated,  as  shown  in 
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Figure  14.  Calculated  current  density 
is  plotted  as  a  function  of  time  for 
diodes  doped  to  1  *  1015  cm2.  Diode 
widths  in  micrometers  are  given  on 
right  of  figure,  whereas  initial  track 
densities  in  cm"3  are  given  on  left. 
Dashed  curves  are  extrapolated. 
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Figure  15.  Log-log  plot  of  time  and 
current  ratios  (defined  in  text)  as 
a  function  of  diode  widths.  Param¬ 
eter  is  initial  track  density  in 

cm"  3. 


figure  15(b).  These  extrapolation  expressions  were  used  to  construct  the 
current-density-versus-time  curves  for  an  initial  track  density  of 
1  x  1018  cm-3  and  diodes  widths  of  1  and  2  pm,  as  shown  by  the  dashed  curves 
of  figure  14.  The  dashed  curves  were  replotted  on  a  linear  scale,  and  the 
total  charge  collected  was  obtained  graphically.  In  each  case  the  total 
charge  was  found  to  be  equal  to  the  initial  charge  in  the  track. 

The  distance  corresponding  to  the  time  ratio,  Tt/TQ  -  4  (an  average 
value  was  chosen),  is  plotted  in  figure  16  on  a  log-log  plot  as  a  function  of 
track  density.  A  straight  line  of  the  slope  -0.75  is  shown.  Likewise,  a 
similar  plot  for  a  fixed  current  ratio  has  the  same  slope.  These  extra¬ 
polations  were  found  to  be  useful  in  estimating  collection  times  for 


calculations  made  with  track  densities  higher  than  1  *  10 


1 8  cm"3 


3 . 2  Gallium  Arsenide 

The  material  parameters  of  gallium  arsenide  are  not  as  well  known  as 
those  of  silicon.  The  ionization  coefficients  were  chosen  to  give  a 
reasonable  fit  to  the  data  of  Bulman  et  al  [9].  Two  electron-velocity-versus- 
field  curves  were  used.  The  parameters  are  listed  in  table  1.  The  hole 
mobility  was  chosen  as  400  cnr/V-s,  and  the  saturation  velocity  of 
1  x  io?  cm/s  attained  at  a  field  of  5  »  101*  V/cm.  The  other  material  param¬ 
eters  were  taken  from  Sze  [6], 


Figure  16.  Distance  required  for 
T  -  4Tq  and  for  J  -  Jq/40  is 
plotted  against  track  density. 


TABLE  1.  ELECTRON  VELOCITY  VERSUS  FIELD  PARAMETERS 
FOR  GaAs  USED  IN  SIMULATIONS  1  AND  2 


Parameter 

Simulation  1 

Simulation  2 

Units 

Low  field  mobility 

7000 

8500 

cm2/V-s 

Maximum  velocity 

1.9 

X 

107 

2.1 

X 

107 

cm/ s 

Field  at  maximum  velocity 

5.8 

X 

103 

5.0 

X 

103 

V/cm 

Minimum  velocity 

6.8 

X 

106 

6.0 

X 

106 

cm/ s 

Field  at  minimum  velocity 

5.0 

X 

104 

1.3 

X 

IQ8 

V/cm 

Saturated  velocity 

1 

X 

107 

1 

X 

107 

cm/s 

The  avalanche  breakdown  voltage  for  a  20-vim  GaAs  diode  doped 
to  1  *  10^5  cm-3  was  found  to  be  295  V,  compared  to  250  V  in  silicon.  Calcu¬ 
lations  have  been  made  for  only  two  track  densities,  1  »  lO1^  cm- 3  and 
1  x  1018  cm~3.  The  collection  currents  as  a  function  of  time  were  quite 
similar  to  those  for  silicon  in  each  case.  These  current  curves  were  also 
nearly  the  same  for  the  two  electron-velocity-versus-f ield  curves,  although 
the  electron  distributions  at  fixed  times  were  greatly  changed.  Figure  17(a) 
shows  the  distribution  of  carriers  at  selected  times  for  a  20-ym  P-type  GaAs 
diode  doped  to  1  x  1018  cm“3.  The  track  density  is  1  *  1014  cm"3,  the  voltage 
is  260  V,  and  the  ionization  coefficients  were  set  to  zero.  The  electron  and 
hole  velocities  are  plotted  as  a  function  of  distance  in  figure  17(b). 


For  the  P-type  diode  the  field  is  low,  and  the  electron  velocity  is 
high  on  the  cathode  (left)  side  of  the  figure.  The  higher  velocity  of  the 
trailing  electrons  causes  the  peaking  of  the  electron  distribution  at  the 
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Figure  17.  (a)  Carrier  distribution  at  selected  times  in  picoseconds  for 
P-type  GaAs  diode  doped  to  1  x  101®  cm"3.  Electron  distributions  are  shown 
as  solid  curves  while  hole  distributions  are  dashed,  (b)  Electron  (solid) 
and  hole  (dashed)  velocities  are  given  as  a  function  of  position. 

trailing  edge.  The  velocity  of  the  trailing  edge  of  the  electron  distribution 
(specifically  the  5  *  101®  cm"  3  location)  is  found  to  be  6  x  io®  cm/s,  the 
input  valley  velocity.  The  holes  accumulate  near  the- cathode  because  their 
velocity  is  low  there.  For  the  complementary  N-type  diode,  the  high  electron 
velocity  region  is  near  the  anode,  and  a  depletion  of  electrons  develops 
there.  Despite  the  marked  differences  in  the  carrier  distributions  between 
the  P-  and  N-type  diodes,  the  maximum  difference  in  current  density  at  any 
time  is  3  percent.  For  the  2  *  10^®  cm"®  track  density,  the  most  noticeable 
effect  of  the  negative-velocity-versus-f leld  characteristic  is  that  the  field 
is  essentially  forbidden  to  be  in  the  range  for  a  negative  velocity  charac¬ 
teristic,  with  abrupt  transitions  in  both  time  and  space.  Again,  there  is 
very  little  effect  on  the  current  density  at  any  time. 


The  recombination  lifetime  in  GaAs  is  not  known  even  in  order  of  mag¬ 
nitude.  Sze  [6]  gives  a  lifetime,  tp,  of  -10"®  s.  In  the  computer  program 
DIODE  the  effective  lifetime  is  inversely  proportional  to  the  carrier 
density.  Calculations  have  been  made  in  GaAs  with  infinite  lifetime,  with 
tR  “  s  and  with  tp  *  2  *  10"9  s  for  N  -  1  *  10^®  cm"®.  For  the  track 
density  1  *  101®  cm"®  calculations,  the  infinite  lifetime  calculations  showed 
collection  currents  closely  agreeing  with  those  of  silicon.  For  the 
t  -  2  *  10"9  s  calculations,  the  total  charge  collected  and  the  time  for  that 
collection  were  reduced  by  approximately  one-half.  This  recombination  effect 
is  larger  than  that  observed  experimentally,  but  avalanching  could  reduce  the 
charge  loss  due  to  recombination.  The  calculations  indicate  that  the  chief 
practical  difference  in  charge  collection  between  silicon  and  GaAs  is  that  the 
latter  has  a  higher  recombination  rate,  reducing  the  collected  charge. 


4.  DISCUSSION 


The  most  Important  question  to  be  discussed  Is  why  our  calculations  show 
avalanche  multiplication  to  be  Important,  whereas  other  calculations  [10]  do 
not.  One  consideration  Is  the  grid  spacing  for  calculations.  In  figure  7  it 
may  be  seen  that  at  the  time  of  the  peak  electric  field  the  entire  high  field 
region  Is  about  0.05-pm  wide.  Certainly  the  grid  spacing  must  be  small  com¬ 
pared  to  this  distance.  Kreskovsky  and  Grubln  [10]  used  28  unequally  spaced 
grid  points  for  their  three-dimensional  calculations  for  a  5-ym  silicon 
diode.  Their  densest  spacing  seems  to  be  about  0.05  pm.  The  maximum  fields 
that  Kreskovsky  and  Grubln  [10]  calculated  were  below  the  avalanche  threshold, 
and  they  did  not  make  avalanche  calculations. 

The  DIODE  computer  program  uses  explicit  difference  equations,  and  for  com¬ 
putational  stability  the  time  step  must  be  less  than  the  dielectric  relaxation 
time,  tD  »  c/a,  where  e  is  the  dielectric  constant  and  o  is  the  conduc¬ 
tivity.  This  Is  actually  a  physical  constraint  for  accurate  space  charge 
calculations.  Further,  the  grid  spacing  must  be  small  compared  to  the 
van  Roosbroeck  [11]  polarization  distance,  X_  -  (v_  +  vn)tp,  where  vp  and  vn 
are  the  hole  and  electron  velocities,  respectively.  When  a  normal ^distri¬ 
bution  of  minority  carriers  is  injected  into  a  semiconductor,  majority 
carriers  largely  neutralize  the  space  charge  in  small  multiples  of  tD,  but  the 
majority  carrier  distribution  is  dragged  along  behind  the  minority  carrier 
distribution  at  a  distance  Xp. 

For  a  track  density  of  1  *  101®  cm"3  in  silicon,  tp  -  3  *  10"1®  s  and  X  - 
6  *  10"®  cm  for  saturation  velocities  of  1  *  10^  cm/s.  The  calculationsPin 
figure  6  were  actually  made  with  a 
timestep  of  5  *  10"15  s  and  a  grid 
distance  of  2  *  10" 7  cm.  A  partial, 
but  recoverable,  instability  is  re¬ 
sponsible  for  the  fluctuations  noted 
in  figure  6.  If  either  the  time  or 
distance  step  is  doubled,  irrecover¬ 
able  instability  results.  A  plot 
comparing  successful  and  failed  runs 
is  shown  in  figure  18.  The  time 
step  and  the  distance  step  used  in 
each  run  are  plotted  as  individual 
points  against  the  track  density  in 
the  figure.  The  straight  lines  are 
the  pertinent  dielectric  relaxation 
times  and  polarization  distance  as  a 
function  of  carrier  density.  Calcul¬ 
ations  made  at  lower  track  densities 
have  shown  that  the  temporal  current 
decay  is  unaffected  by  the  recover¬ 
able  instability.  Using  these  time 
and  distance  steps  in  a  two-  or 
three-dimensional  (2D  or  3D)  calcu¬ 
lation  will  be  quite  costly. 


Figure  18.  Dielectric  relaxation  time 
(solid  line)  and  polarization  distance 
(dashed  line)  are  shown  as  a  function 
of  carrier  density.  Individual  points 
are  time  and  distance  steps  which  re¬ 
sulted  in  successful  or  failed  computer 
runs . 


Some  indication  of  the  2D  aspect  of  charge  collection  can  be  deduced  from 
ID  calculations.  Figure  19(a)  shows  the  field  distribution  calculated  at  42 
ps  for  a  track  density  of  2  *  10^5  cm-3  in  a  20-um  N-type  diode  doped  to 
1  *  1015  cra"3.  The  applied  voltage  is  240  V  and  the  avalanche  constants  were 
set  to  zero.  Electrons  are  depleted  from  the  first  4  ym,  increasing  the  field 
gradient,  and  holes  from  the  last  2  pm,  reversing  the  field  gradient.  The 
latter  is  an  example  of  field  funnel ing.  Figure  19(b)  shows  a  sketch  of 
constant  field  contours  deduced  from  figure  19(a).  The  vertical  axis  (not 
Quantified)  is  the  track  radius  with  the  track  centered.  The  field  contours 
at  the  edges  are  equidistant,  given  by  the  field  gradient  at  t  »  0.  The  field 
at  the  center  of  the  track  is  determined  by  the  42-ps  curve  of  figure  19(a). 
The  full  contour  is  drawn  freehand  from  the  center  to  the  end  points.  Current 
paths  will  follow  the  field  lines,  which  are  perpendicular  to  the  field  con¬ 
tours.  Thus,  the  electrons  will  diverge  from  the  track  in  the  first  3  pm  and 
last  6  um  from  the  cathode,  and  converge  in  between.  For  this  example  the 
maximum  field  at  the  cathode  was  390  kV/cm;  the  maximum  was  reduced  to 
310  kV/cm  for  calculations  including  avalanche  multiplication.  Some  calcu¬ 
lations  at  a  track  density  of  1  «  101®  cm"3  show  higher  fields  with  ava¬ 
lanching  than  those  without. 

Diffusion  in  one  and  two  dimensions  can  also  be  used  to  obtain  some 
estimate  of  the  error  of  ID  as  compared  with  2D  calculations.  The  peak 
density  of  2D  diffusion  from  a  line  source  decays  at  the  same  rate  as  the  ID 
diffusion  from  a  plane  source  with  double  the  diffusion  constant.  From  this 
analogy  and  from  the  consideration  of  the  previous  paragraph,  an  estimate  is 
made  that  the  2D  a  calculation  would  have  a  maximum  field  of  about  one-half 
that  calculated  for  the  ID  calculation.  In  this  context  our  ID  calculation 
would  have  about  the  same  avalanche  threshold  and  charge  multiplier  as  would 
be  obtained  in  a  2D  calculation  with  double  the  track  density.  When  it  is 
considered  that  the  track  cross-sectional  area  may  be  unknown  by  an  order  of 
magnitude,  the  ID  calculation  seems  tolerable. 
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Figure  19.  Field  distributions  for 
a  20-pm  N-type  diode  doped  to 
1  x  1015  cm-3  and  an  initial  track 
density  of  2  x  lO^  cm-3;  (a)  Solid 
line  is  field  at  42  ps;  dashed  line 
is  for  t  »  0  and  (b)  Two-dimensional 
field  contour  lines  deduced  from  (a). 


Figure  3  shows  close  agreement  between  Vmln,  the  minimum  voltage  for  ava¬ 
lanche  multiplication  as  a  function  of  track  density,  and  VB0  ,  the  avalanche 
breakdown  of  one-sided  diodes  as  a  function  of  doping  density,  for  carrier 
densities  above  lO1^  cm"3.  The  deviation  below  this  density  is  due  to  the 
device  width  of  20  pm  for  the  track  calculation.  This  agreement  may  be  under¬ 
stood  from  two  considerations.  First,  in  the  one-sided  diodes,  avalanche 
multiplication  is  noted  at  about  one-half  of  VBQ.  Second,  the  high  field 
region  near  the  cathode  is  almost  identical  in  the  two  calculations,  since 
depletion  is  essentially  complete  in  each  case.  However,  in  the  track  cal¬ 
culation,  there  is  also  a  high  field  region  near  the  anode  due  to  hole 
depletion  (see  fig.  7).  Therefore,  the  total  voltage  is  nearly  the  same  in 
each  case. 

It  was  shown  in  figure  7  that  the  maximun  field  was  obtained  in  less  than 
1  ps.  This  time  is  comparable  to  the  time  for  the  ionizing  particle  to  com¬ 
pletely  form  the  track.  Therefore,  an  accurate  calculation  should  include 
the  formation  of  the  track  with  the  velocity  of  the  ionizing  particle. 

The  mobilities  used  in  the  calculations  were  characteristic  of  the 
1015  cm“3  doping  level  rather  than  the  usually  higher  track  density.  This 
assumption  has  the  effect  of  increasing  the  current  and  thus  reducing  the  time 
scale  of  the  current  pulse.  However,  carriers  were  at  a  saturated  velocity 
for  the  most  part,  and  saturated  velocities  are  independent  of  carrier 
density. 

The  closeness  of  the  agreement  between  the  ID  calculation  result  and  the 
measured  data  3hown  in  figure  5  was  unexpected.  Various  compensating  errors 
may  account  for  this.  Campbell  et  al  [2]  estimate  their  diode  width  as  1  pm, 
but  their  breakdown  voltage  was  about  7  V.  Our  calculations  show  that  a  1-pm 
diode  should  have  a  breakdown  voltage  at  about  35  V.  Although  calculated 
breakdown  voltages  often  exceed  measured  ones  for  various  reasons,  one 
suspects  that  their  diode  is  no  more  than  0.5  pm.  If  this  were  so,  their 
charge  collection  efficiency  would  be  -1.0  instead  of  -0.5,  as  they  give  for 
zero  applied  voltage.  This  was  assumed  in  our  figure  5.  Our  calculations 
assume  that  the  ionized  track  was  fully  formed  at  zero  time,  but  as  seen  in 
the  previous  paragraph,  the  maximum  field  was  obtained  in  a  time  comparable  to 
the  transit  time  of  the  ionizing  particle.  Therefore,  the  high  field  at  the 
anode  due  to  depletion  of  holes  would  be  delayed,  and  more  voltage  would  be 
available  for  avalanche  multiplication  near  the  cathode.  This  possible  factor 
of  two  field  enhancement  would  negate  the  postulated  factor  of  two  decrease 
due  to  the  ID  calculations. 

5.  CONCLUSIONS 

One-dimensional  calculations  show  that  avalanche  multiplication  must  be 
considered  in  the  charge  collection  in  single  ionizing  events.  The  calculated 
increase  in  charge  multiplication  as  a  function  of  voltage  is  shown  to  be  in 
good  agreement  with  the  measurements  of  Campbell  et  al  [2].  Charge  multi¬ 
plication  becomes  more  important  as  the  epitaxial  layer  decreases  in 
thickness.  Since  one  hardening  technique  is  to  reduce  field  funneling  by 
using  thinner  epitaxial  layers,  one  must  beware  of  increasing  the  charge  mul- 
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tiplicatlon  by  a  large  factor  due  to  avalanching  while  avoiding  a  limited 
multiplication  due  to  Tunneling. 
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